Dynamics and energetics of bubble growth in magmas: Analytical formulation and numerical modeling

We have developed a model of diffusive and decompressive growth of a bubble in a finite region of melt which accounts for the energetics of volatile degassing and melt deformation as well as the interactions between magmatic system parameters such as viscosity, volatile concentration, and diffusivity. On the basis of our formulation we constructed a numerical model of bubble growth in volcanic systems. We conducted a parametric study in which a saturated magma is instantaneously decompressed to one bar and the sensitivity of the system to variations in various parameters is examined. Variations of each of seven parameters over practical ranges of magmatic conditions can change bubble growth rates by 2–4 orders of magnitude. Our numerical formulation allows determination of the relative importance of each parameter controlling bubble growth for a given or evolving set of magmatic conditions. An analysis of the modeling results reveals that the commonly invoked parabolic law for bubble growth dynamics R∼t1/2 is not applicable to magma degassing at low pressures or high water oversaturation but that a logarithmic relationship R∼log(t) is more appropriate during active bubble growth under certain conditions. A second aspect of our study involved a constant decompression bubble growth model in which an initially saturated magma was subjected to a constant rate of decompression. Model results for degassing of initially water‐saturated rhyolitic magma with a constant decompression rate show that oversaturation at the vent depends on the initial depth of magma ascent. On the basis of decompression history, explosive eruptions of silicic magmas are expected for magmas rising from chambers deeper than 2 km for ascent rates \u3e1–5 m s−1

Dynamics of coupled diffusive and decompressive bubble growth in magmatic systems

Bubble growth in an ascending parcel of magma is controlled both by diffusion of oversaturated volatiles and decompression as the magma rises. We have developed a numerical model which explores the processes involved in water exsolution from basaltic and rhyolitic melts rising at a constant rate from magma chamber depths of 4 and 1 km. While the model does not attempt to simulate natural eruptions, it sheds light on the processes which control eruptive behavior under various conditions. Ascent rates are defined such that a constant rate of decompression dP/dt is maintained. A variety of initial ascent rates are considered in the model, from 1 m/s to 100 m/s for basalts, and from a few centimeters per second to 10 m/s for rhyolite, at the base of the conduit. The model results indicate that for any reasonable ascent rate, basaltic melt degasses at a rate sufficient to keep the dissolved volatile concentration at equilibrium with the decreasing ambient pressure. Rhyolitic melt reaches the surface at equilibrium if its ascent rate is less than 1 m/s, but it can erupt with high oversaturation at greater ascent rates. The latter may lead to explosive eruptions. If the ascent rate of rhyolite is 10 m/s or more, then melt barely degasses at all in the conduit and erupts with the highest oversaturation possible. For the case of slow magma rise, bubble growth is limited by decompression. For the case of rapid magma rise, bubble growth is limited by diffusion. The results of our simple model do not accurately simulate natural volcanic eruptions, but suggest that subsequent, more complex models may be able to simulate eruptions using the insights regarding diffusive and decompressive bubble growth processes explored in this study. Numerical modeling of volcanic degassing may eventually lead to better prediction of eruption timing, energetics and hazards of active volcanoes

Reply [to “Comment on “Dynamics and energetics of bubble growth in magmas: Analytical formulation and numerical modeling” by A. A. Proussevitch and D. L. Sahagian”]

We have developed a model of diffusive and decompressive growth of a bubble in a finite region of melt which accounts for the energetics of volatile degassing and melt deformation as well as the interactions between magmatic system parameters such as viscosity, volatile concentration, and diffusivity. On the basis of our formulation we constructed a numerical model of bubble growth in volcanic systems. We conducted a parametric study in which a saturated magma is instantaneously decompressed to one bar and the sensitivity of the system to variations in various parameters is examined. Variations of each of seven parameters over practical ranges of magmatic conditions can change bubble growth rates by 2–4 orders of magnitude. Our numerical formulation allows determination of the relative importance of each parameter controlling bubble growth for a given or evolving set of magmatic conditions. An analysis of the modeling results reveals that the commonly invoked parabolic law for bubble growth dynamics R∼t1/2 is not applicable to magma degassing at low pressures or high water oversaturation but that a logarithmic relationship R∼log(t) is more appropriate during active bubble growth under certain conditions. A second aspect of our study involved a constant decompression bubble growth model in which an initially saturated magma was subjected to a constant rate of decompression. Model results for degassing of initially water‐saturated rhyolitic magma with a constant decompression rate show that oversaturation at the vent depends on the initial depth of magma ascent. On the basis of decompression history, explosive eruptions of silicic magmas are expected for magmas rising from chambers deeper than 2 km for ascent rates \u3e1–5 m s−1

Environmental science applications with Rapid Integrated Mapping and analysis System (RIMS)

The Rapid Integrated Mapping and analysis System (RIMS) has been developed at the University of New Hampshire as an online instrument for multidisciplinary data visualization, analysis and manipulation with a focus on hydrological applications. Recently it was enriched with data and tools to allow more sophisticated analysis of interdisciplinary data. Three different examples of specific scientific applications with RIMS are demonstrated and discussed. Analysis of historical changes in major components of the Eurasian pan-Arctic water budget is based on historical discharge data, gridded observational meteorological fields, and remote sensing data for sea ice area. Express analysis of the extremely hot and dry summer of 2010 across European Russia is performed using a combination of near-real time and historical data to evaluate the intensity and spatial distribution of this event and its socioeconomic impacts. Integrative analysis of hydrological, water management, and population data for Central Asia over the last 30 years provides an assessment of regional water security due to changes in climate, water use and demography. The presented case studies demonstrate the capabilities of RIMS as a powerful instrument for hydrological and coupled human-natural systems research

Development of Distributed Research Center for analysis of regional climatic and environmental changes

We present an approach and first results of a collaborative project being carried out by a joint team of researchers from the Institute of Monitoring of Climatic and Ecological Systems, Russia and Earth Systems Research Center UNH, USA. Its main objective is development of a hardware and software platform prototype of a Distributed Research Center (DRC) for monitoring and projecting of regional climatic and environmental changes in the Northern extratropical areas. The DRC should provide the specialists working in climate related sciences and decision-makers with accurate and detailed climatic characteristics for the selected area and reliable and affordable tools for their in-depth statistical analysis and studies of the effects of climate change. Within the framework of the project, new approaches to cloud processing and analysis of large geospatial datasets (big geospatial data) inherent to climate change studies are developed and deployed on technical platforms of both institutions. We discuss here the state of the art in this domain, describe web based information-computational systems developed by the partners, justify the methods chosen to reach the project goal, and briefly list the results obtained so far

Analysis of Vesicular Basalts and Lava Emplacement Processes for Application as a Paleobarometer/Paleoaltimeter

We have developed a method for determining paleoelevations of highland areas on the basis of the vesicularity of lava flows. Vesicular lavas preserve a record of paleopressure at the time and place of emplacement because the difference in internal pressure in bubbles at the base and top of a lava flow depends on atmospheric pressure and lava flow thickness. At the top of the flow, the pressure is simply atmospheric pressure, while at the base, there is an additional contribution of hydrostatic lava overburden. Thus the modal size of the vesicle (bubble) population is larger at the top than at the bottom. This leads directly to paleoatmospheric pressure because the thickness of the flow can easily be measured in the field, and the vesicle sizes can now be accurately measured in the lab. Because our recently developed technique measures paleoatmospheric pressure, it is not subject to uncertainties stemming from the use of climate‐sensitive proxies, although like all measurements, it has its own sources of potential error. Because measurement of flow thickness presupposes no inflation or deflation of the flow after the size distribution at the top and bottom is “frozen in,” it is essential to identify preserved flows in the field that show clear signs of simple emplacement and solidification. This can be determined by the bulk vesicularity and size distribution as a function of stratigraphic position within the flow. By examining the stratigraphic variability of vesicularity, we can thus reconstruct emplacement processes. It is critical to be able to accurately measure the size distribution in collected samples from the tops and bottoms of flows because our method is based on the modal size of the vesicle population. Previous studies have used laborious and inefficient methods that did not allow for practical analysis of a large number of samples. Our recently developed analytical techniques involving high‐resolution x‐ray computed tomography (HRXCT) allow us to analyze the large number of samples required for reliable interpretations. Based on our ability to measure vesicle size to within 1.7% (by volume), a factor analysis of the sensitivity of the technique to atmospheric pressure provides an elevation to within about ±400 m. If we assume sea level pressure and lapse rate have not changed significantly in Cenozoic time, then the difference between the paleoelevation “preserved” in the lavas and their present elevation reflects the amount of uplift or subsidence. Lava can be well dated, and therefore a suite of samples of various ages will constrain the timing of epeirogenic activity independent of climate, erosion rates, or any other environmental factors. We have tested our technique on basalts emplaced at known elevations at the base, flanks, and summit of Mauna Loa. The results of the analysis accurately reconstruct actual elevations, demonstrating the applicability of the technique. The tool we have developed can subsequently be applied to problematic areas such as the Colorado and Tibetan Plateaus to determine the history of uplift

The use and re-use of unsustainable groundwater for irrigation: A global budget

Depletion of groundwater aquifers across the globe has become a significant concern, as groundwater is an important and often unsustainable source of irrigation water. Simultaneously, the field of water resource management has seen a lively debate over the concepts and metrics used to assess the downstream re-use of agricultural runoff, with most studies focusing on surface water balances. Here, we bring these two lines of research together, recognizing that depletion of aquifers leads to large amounts of groundwater entering surface water storages and flows by way of agricultural runoff. While it is clear that groundwater users will be impacted by reductions in groundwater availability, there is a major gap in our understanding of potential impacts downstream of groundwater pumping locations. We find that the volume of unsustainable groundwater that is re-used for irrigation following runoff from agricultural systems is nearly as large as the volume initially extracted from reservoirs for irrigation. Basins in which the volume of irrigation water re-used is equal to or greater than the volume of water initially used (which is possible due to multiple re-use of the same water) contain 33 million hectares of irrigated land and are home to 1.3 billion people. Some studies have called for increasing irrigation efficiency as a solution to water shortages. We find that with 100% irrigation efficiency, global demand for unsustainable groundwater is reduced by 52%, but not eliminated. In many basins, increased irrigation efficiency leads to significantly decreased river low flows; increasing irrigation efficiency to 70% globally decreases total surface water supplies by backsim600 km3 yr−1. These findings illustrate that estimates of aquifer depletion alone underestimate the importance of unsustainable groundwater to sustaining surface water systems and irrigated agriculture